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One of the promising systems to realize quantum computation is a hybrid system where a superconducting 
flux qubit plays a role of a quantum processor and the NV center ensemble is used as a quantum memory. We 
have theoretically and experimentally studied the effect of magnetic fields on this hybrid system, and found that 
the lifetime of the vacuum Rabi oscillation is improved by applying a few mT magnetic field to the NV center 
ensemble. Here, we construct a theoretical model to reproduce the vacuum Rabi oscillations with/without 
magnetic fields applied to the NV centers, and we determine the reason why magnetic fields can affect the 
coherent properties of the NV center ensemble. From our theoretical analysis, we quantitatively show that the 
magnetic fields actually suppress the inhomogeneous broadening from the strain in the NV~ centers. 


Hybridization is a promising approach for quantum compu¬ 
tation Cli. Each system has characteristic with its own ad¬ 
vantages and disadvantages. To couple different systems, we 
hope to pick up the advantage which each system has. One 
of the candidates for such hybrid systems is a superconduct¬ 
ing circuit such as a superconducting flux qubit (FQ) and an 
electron spin ensemble such as nitrogen-vacancy (NV“) cen¬ 
ters 1131 4201] . as described in the Fig. [T] High controllability 
of superconducting FQs has already been achieved with exist¬ 
ing technology 112111 . Reliable gate operations have been al¬ 
ready demonstrated 11221] . Quantum non-demolition measure¬ 
ments can be performed by Josephson bifurcation amplifier 
Ein . However, despite significant effort, the coherence time 
of the FQ is of the order of 10 ps 23]. On the other hand, 
the NV“ center has a long coherence time |24-^. With 
dynamical decoupling, the coherence time of electron of the 
NV“ center is 0.6 s ]|29[] that is much longer than the FQ. So, 
coupling the FQ with the NV“ centers is a promisiim way to 
obtain both controllability and long coherence time 1314201]. 

It is often useful to transfer the state between the FQ 
and NV“ centers for quantum information processing. We 
keep quantum states in the quantum memory (NV“ centers) 
when gate operations are not required. On the other hand, 
to perform gate operations, we need to transfer the quantum 
states from the quantum memory to the quantum processor 
(FQ), which can be realized by using vacuum Rabi oscillation 
(VRO). However, the error rate of the state transfer in the cur¬ 
rent technology is an order of ten percent [Toi [b^, which 
is too large to perform quantum computation 011 1321] . The 
noise mainly comes from the inhomogeneous broadening of 
the NV“ centers. Therefore, it is crucial to suppress the deco¬ 
herence of the NV“ centers for computational tasks. 

Improving the coherence time of the NV“ center ensem¬ 
ble also has a fundamental importance in the area of quantum 
metrology Sin, quantum walk S, and quantum simu¬ 
lation 1351]. In these applications, the efficiency strongly de¬ 
pends on the coherence time of the ensemble of NV“ centers. 
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FIG. 1: Illustration of the hybrid system composed of a supercon¬ 
ducting flux qubit and an ensemble of NV“ centers. The flux qubit 
has four Josephson junctions that forms a two-level system. There 
are two control lines for the flux qubit. We use one of the control 
lines to change the energy bias e, and use the other line to change the 
energy gap A. Diamond crystal is glued on top of a flux qubit, and 
this diamond contains NV~ centers. The electrons spins trapped in 
the NV~ center form a three-level system, and so we have a V-type 
energy level structure for the NV“ center. 


So it is essential in these areas to And a way to improve the 
coherence time of the NV“ centers. 

A cavity protection ][I8, 36, 37] is a promising way to im¬ 
prove the coherence time of the NV“ center ensemble. If the 
coupling strength between the ensemble of NV“ centers and 
a superconducting flux qubit (or microwave cavity) is larger 
than the inhomogeneous width of the NV“ centers, the col¬ 
lective mode of the NV“ centers could be well decoupled 
from the other sub-radiant states of the NV“ so that inho¬ 
mogeneous broadening would be suppressed. However, it has 
been shown that, if the spectral density of the inhomogeneous 
broadening is described by a Forentzian, the noise cannot be 
suppressed by the cavity protection OlfllHOT . Moreover, for the 
applications of quantum memory and quantum held sensing. 
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it is sometime necessary to turn off the interaction between 
the NV centers and superconducting circuit. So it is better to 
have an alternative scheme that will work even for such cases. 


In this paper, we report an improvement of the lifetime of 
the VRO by applying an in-plain magnetic field to this hybrid 
system. We have observed VRO with/without the magnetic 
field, and the lifetime of the VRO with magnetic field is nearly 
twice than without magnetic field. We have constructed a the¬ 
oretical model to reproduce these results, and have found that 
the magnetic field suppress the inhomogeneous strain effect. 

Let us now describe our experimental set-up. Our system 
consists of a gap-tunable FQ ^ on which a diamond 
crystal with NV“ density of approximately 5 x 10^^ cm“^ 
is bonded d. The NV ensemble is created by the ion 
implantation and annealing in in vacuum 1 13 3l- The dis¬ 
tance between the FQ and the surface of the diamond crystal 
is less than 1/rm. We can apply external magnetic field of 
2.6 mT along the [100] crystalline axis lll6l] . The gap-tunable 
FQ is fabricated with a superconducting loop containing four 
Josephson junctions. To control the FQ, a microwave line is 
fabricated around the FQ. The FQ is designed to couple with 
the SQUID structure via magnetic fields. The probability of 
the excited state of the FQ is measured by the SQUID. 

To observe the VRO, we perform the following experiment. 
First, we excited the FQ by applying microwave pulse where 
the FQ is decoupled from the NV“ centers by the detuning. 
Second, we brought the FQ into the resonance of the NV— 
centers by applying a magnetic flux. Finally, after a time t, we 
can measure the excited probability of the FQ via the SQUID. 
The measurements were done in a dilution refrigerator at a 
temperature below 50 mK. 

We describe our system by the Hamiltonian ^ , 13, 18. 191 
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where ax,y,z denotes the Pauli matrix for FQ with cTx whose 
eigenstates corresponds to two persistent-current states. Also, 
we define (T± = ^{^x ± The electron spin 1 opera¬ 

tors of the NV“ center are represented by Sx^y^z- T^flux de¬ 
notes the Hamiltonian of the FQ where A denotes the en¬ 
ergy gap and e denotes the magnetic energy bias. iFens rep¬ 
resents the ensemble composed of N individual NV“ cen¬ 
ters where Z)/27r ~ 2.878GHz denotes a zero field splitting, 

Ek = \J{eI'^Y + (-£' 2 ^^)^ denotes a strain induced splitting, 
geMsBNv ■ S denotes a Zeeman splitting, and Bnv denotes 


a magnetic field with pe/iB/27r = 28MHz-mT~^. A quan¬ 
tization axis (z axis) to be the direction from the vacancy to 
the nitrogen is set by the zero field splitting of the NV“ cen¬ 
ter. For a small magnetic field D ^ pg/rsIBNvl, the x and 
y component of the magnetic field is insignificant to change 
quantized axis of the NV“ center, and so we include only 
the effect of z axis of the field. We consider three relevant 
type of the magnetic field Bnv: an in-plane external magnetic 
field Bex, an inhomogeneous magnetic field due to PI centers 
Binh, and a hyperfine field from the nitrogen nuclear spins Bhf. 
The term Elint denotes the magnetic coupling between the FQ 
and the NV“ centers where represented the magnetic 
field induced by persistent current of the FQ. Since collec¬ 
tive enhancement of the coupling strength between the NV“ 
and FQ is not available along the z axis of the NV“ center 
11311, we can ignore the coupling with Sz- We can write Hint 
as Hint = ^Y^k=i9^x ' {Sx,k cos (j)k - Sy,k sin (j)k) where 


g = denotes a Zeeman splitting of the NV spin 

due to FQ magnetic field in x-y plane and tpk denotes 

the angle of the field in the plane. 

For an ensemble of NV^ center with only a few excitations 
in it at most, we can use the Holstein-Primakoff approxima¬ 
tion to treat NV“ spins as an ensemble of harmonic oscillators 
0. We define creation (destruction) operators of a bright 
state and a dark state of the NV center lfl3l by bj, = |A+)fc(0| 
and 4 = |^-)fe(0|, {bk = |0)fc(A+| and 4 = |0)fe(A_|) 
where \A+)k = ^(|l)fc + |-l)fe) + ^(|l)fc-|-l)fe), 
\A.)k = ^(|l)fe + |-l)fc) + ^(|l)fe-|-l)fe). Here, 
|A+)fe(0| denotes the state of the NV“ center to be directly 
coupled with the FQ while |A_)fe(0| has no direct coupling 
with the FQ. Since we consider only one or zero excitation in 
a total system, we can also replace a ladder operator of the FQ 
with a creation operator of a harmonic oscillator as it+ —c^. 

Moving to a rotating frame with angular frequency cu 
defined hy U = ^nd making the rotating 

wave approximation, we obtain the simplified Hamiltonian 

H ~ + {^ll"’b\hu + hull'’dldk + hg'{c^bk + 


cbl) + h{Jk + iJ'k)b\^k + HJk - where Wc = 

\/e'^ + A^, LuH'’ ~ Bfe - w - Ey^k, ~ Dk - UJ + Ey^k, 

Jk = geEsBl'l 4 = Ex,k, g' = g^j'Je^ + A^, Ex,k = 

Ek cos 24 , and Ey^k = Ek sin 24 . 

Now the dynamics of this hybrid system can be investigated 
using the Heisenberg equations of motions. We write Heisen¬ 
berg equations of motion as 


d „ 

—c = -iWcC 
dt 


N 

gksin^- bk) - FcC 


(6) 


(X ^ xy A. ^ 

— bk = -iujl ’bk - iJkdk + Jkdk - igsin^k ■ c- Fb^fc (7) 

^dk =-iujf'’dk - iJkh - J'kh-^ddk (8) 

where Fc, F;,, and F^ denote the decay rate of c, b, and d, re- 
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spectively. We numerically solve these equations with an ini¬ 
tial state of |'0(t = 0)) = c^jvac), and plot the renormalized 
excitation probability of the FQ (which corresponds a switch¬ 
ing probability of SQUID) in the Fig. |2]and Fig. |3 



Time (ns) 

FIG. 2: Vacuum Rabi oscillations between the FQ and NV~ cen¬ 
ters without an applied external magnetic field. Red dots show the 
experimental results while the blue line shows the results of a nu¬ 
merical model where we use N = 1200, e = 0, rc/27r = 0.3MHz, 
r6/27r = r<i/27r = 0.44MHz, JDfc/27r = 0.08MHz (FWHM), 
S{gfj.BBi’''>)/2-K = 3.1MHz(FWHM), SEi^kl^n = 5E2,kl2-K = 
4.4MHz(FWHM), Ahf ^ 2.3MHz, and ^/kg/2n = 13MHz 
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FIG. 3: Vacuum Rabi oscillations between the FQ and NV“ centers 
with an applied external magnetic field of 2.6 mT along the [100] 
crystalline axis. Except the magnetic field, we use the same parame¬ 
ters as those in Fig. |2] 

Now let us detail our simulation technique and the core el¬ 
ements of it. The Lorentzian distributions, which have been 
typically used to describe the inhomogeneous broadening of 
the NV“ centers |@,[l3,[I3>li3]> are assumed for Ei,k, 
and E 2 ,k {k = - ,N) to include the effect of the inho¬ 

mogeneous lattice distortion of the NV“ centers. Next, due to 
the electron spin-half bath in the environment such as PI cen¬ 
ter, randomized magnetic field on the NV center exists, and 
the nitrogen nuclear spin splits the electron-spin energy into 
three level via a hyperfine coupling. To include both these 
two effect, we use a random distribution of the magnetic fields 
with the form of the mixture of three Lorentzian functions that 
are separated with 27r x 2.3 MHz due to the hyperfine inter¬ 
action with nuclear spin iSES. With these assumptions. 


we have reproduced the VRO with/without applied magnetic 
field in the Fig. |2]and Fig. |3] While the FQ can induce the 
transition both |0)fe O | — l)k and |0)fc O |l)fc with zero 
applied magnetic field, the FQ induces only one of them with 
applied magnetic field of a few mT due to the detuning ef¬ 
fect 0113, foi [13,[T^ . This changes the effective coupling 
strength between the FQ and NV centers, and is the cause of 
the different time interval of the oscillations in the Fig. |2]and 
Fig. [ 3 ] With an applied magnetic field of 2.6 mT, the VRO 
can be observed until around 170 ns while we cannot observe 
a clear oscillation beyond 100 ns without an applied magnetic 
field. Thus, it shows the external magnetic field improves the 
lifetime of the VRO. 


We explain how we determine the parameters for the nu¬ 
merical simulations. Fc can be determined by the Ti mea¬ 
surement of the FQ, which was performed independent of the 
VRO experiment |[l3l- Since the frequency shift of E is 50 
times larger than that of D when an electric field is applied 
I 4in . we use 5E ~ 50(5in this paper. It is known that, from 
the spectroscopic measurements, a sharp peak located in the 
middle of the avoided crossing structure was observed in this 
hybrid system 0, and one can determine the r;i(= F^) from 
the width of this sharp peak lfl3] . The time period of the VRO 
let us specify the value of y/iV g. Moreover, as we describe 
later, the envelope of the VRO with (without) magnetic field 
is mainly determined by SB^ and 6E). So, by fitting the 
spectroscopy and VRO with/without magnetic field, one can 
specify necessary parameters for our model. 


We understand why the applied magnetic field actually 
improves the coherence time of the NV~ centers. There 
are two relevant decoherence source for the NV centers, 
inhomogeneous magnetic fields and the strain distribution 
||3 113 O- The magnetic-field noise comes from an ef¬ 
fective Hamiltonian between PI center and NV center as 
Heff = J2k i where flip-flop term is negligi¬ 

ble due to the large energy difference between them. How¬ 
ever, this term commutes with the Zeeman term of the ap¬ 
plied external magnetic field as QefJ-BBexSz, and so exter¬ 
nal magnetic field cannot affect this. On the other hand, the 
Hamiltonian of the inhomogeneous strain does not commute 
with the Zeeman Hamiltonian of the applied external mag¬ 
netic field. Actually, the Hamiltonian of the k th NV“ cen¬ 
ter is written as -F -F ((Jfe -F 

iJ'j^)b\dk + e.c), and the eigenenergies are i?o = 0 and 

E± = D ± Elj^ + E'^ f, -F {geUsBi'"'')'^. If the magnetic 
field is large, we can expand the eigenenergies of the excited 
states as E+ ~ geUsBi^^ ± Sq the effect of the 

variations of E^^k and Ey^k becomes negligible, and this could 
improve the lifetime of the VRO H. 


To confirm this effect, we performed another numerical 
simulation of the VRO with applied magnetic field for the 
strain values of 5E/2tt — 4.4,6.0, 7.6 MHz in the Fig. |4] 
Interestingly, the three VROs shown in the Fig. |4] are almost 
the same. These results clearly show that the lifetime of the 
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VRO is quite robust against the inhomogeneous strains. So we 
conclude that the external magnetic field suppress the effect of 
inhomogeneous strain so that the improvement of the lifetime 
has been observed in our experiment. It is worth mentioning 
that, although possible suppression of the strain distribution 
by the applied magnetic field is mentioned in 14211 . we firstly 
demonstrate such a mechanism by the experiment. 


mogeneous broadening of the strain. This result is a relevant 
step toward to the realization of the long-lived quantum mem¬ 
ory for a superconducting flux qubit. 

We thank R. Amsiiss and H. Nakano for valuable discus¬ 
sions. This work was supported by KAKENHI(S) 25220601 
and in part by Commissioned Research of NICT. 
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FIG. 4: Numerical simulations of vacuum Rabi oscillations with the 
applied magnetic field for several values of the strains. Except the 
values of the strains, we use the same parameters as those in Fig. 


It is known that applying transversal magnetic field can im¬ 
prove the coherence time of the NV“ centers when the strain 
distribution is much larger than the decoherence rate due to 
the randomized magnetic fields Here, the application of 
the transversal magnetic field can suppress the decoherence 
due to the environmental magnetic field while this cannot sup¬ 
press the inhomogeneous broadening of the strain 11431] . On 
the other hand, our scheme to apply horizontal magnetic field 
is complemental to this, because we can suppress the strain 
inhomogeneous broadening. Our scheme has an advantage 
especially when we can decrease the broadening due to the 
magnetic field by using another technique. 

Actually, it is possible to combine our scheme with another 
technique to reduce the inhomogeneous magnetic fields. For 
example, one way to suppress the magnetic noise for the NV“ 
centers is to decrease the PI centers (a nitrogen atom substitut¬ 
ing a carbon atom) by using differently synthesized diamond 
crystals US. PI centers are considered as a electron spin-half 
bath, and they cause randomized magnetic field to decoher the 
NV“ centers SGI. However, reduction of the PI centers 
cannot suppress the noise due to the strain. Therefore, by ap¬ 
plying a magnetic field with a diamond having less PI centers, 
further improvement of the coherence time should be possible. 
This provides us with a sensitive diamond-based field sensor, 
and will be useful for a long-lived quantum memory during 
quantum computation as future applications. 

In conclusion, we have observed an improvement of a life¬ 
time of the vacuum Rabi oscillations between the FQ and NV 
centers by applying an in-plain magnetic field. By reproduc¬ 
ing the experimental result from a theoretical model, we have 
found that the applied magnetic field can suppress the inho- 
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